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Abstract. Population monitoring is a crucial method for conservation projects, especially for the highly endangered clade 
of amphibians which is threatened by habitat loss and emerging infectious diseases such as the chytrid fungus Batracho- 
chytrium salamandrivorans (Bsal). Bsal has led to massive population declines of European fire salamanders in Belgium, 
the Netherlands and parts of Germany and population trends are decreasing. Thus, regular population monitoring is es- 
sential to keep track of population dynamics and detect potential Bsal outbreaks, especially since valid population esti- 
mates for many species and/or populations are scarce. In recent years, photographic mark-recapture studies have gained 
importance enabling researchers to keep track of individuals without the use of harmful marking techniques that might 
influence behaviour or survival. We monitored a European fire salamander population over four years in the Kottenforst 
forest, Germany, by combining a mark-recapture approach with a new photo-recognition software, the Amphibian and 
Reptile Wildbook. We investigated potential differences of two ecotypes, i.e., larval populations bred in ponds and streams. 
Furthermore, we compared the observed number of larvae and the estimated population sizes between two consecutive 
years. There was a year effect on the observed number of larvae, the percentage of injured larva and estimated population 
size. The habitat type affected apparent survival rates. There was also a habitat effect on the percentage of injured larvae, 
but this was only existent in some years. The mean larval size, water temperature and the superpopulation size were not 
affected by any of the aforementioned factors. Although the population was free from Bsal, the high variation in the popu- 
lation estimates emphasize the need for a regular and standardised monitoring to assess the current population status and 
detect early population declines that might otherwise remain overlooked. 
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photographic mark-recapture. 


Introduction 


The lack of reliable population estimates for most of the 
species is one of the major problems in conservation bio- 
logy. We are in the midst of a global biodiversity crisis and 
population trends are decreasing globally (reviewed in 
SINGH 2002). The main threats of biodiversity are changes 
in land and sea use, species overexploitation, pollution, cli- 
mate change, invasive species, and new diseases (SALAF- 
Sky et al. 2008, WWE 2020). Almost one third of species 
worldwide are threatened with extinction (IUCN 2021) 
and studies indicate that we are close to a sixth mass ex- 
tinction (e.g., WAKE & VREDENBURG 2008, BARNOSKY et al. 
2011, MCCALLUM 2015). 

Amphibians, with 41% endangered species (IUCN 
2021), are the most endangered group (STUART et al. 2004), 
but have gained far less attention than other taxa (reviewed 
in BRITO 2008). The causes of amphibian population de- 
clines are complex, spatially different (SopHr et al. 2008, 
GRANT et al. 2016) and include multiple factors such as cli- 


matic or environmental conditions and the spread of new 
diseases (e.g., COLLINS 2010, MANENTI et al. 2022). 

Thus, monitoring efforts have increased over the last 
decades (reviewed in Moussy et al. 2022) for purposes such 
as status assessment, measuring the effectiveness of con- 
servation approaches or the estimation of population sizes 
(e.g., DONNELLY & GUYER 1994, STEM et al. 2005). Popula- 
tion size estimation for open populations requires the indi- 
vidual identification of animals, which, in past studies, was 
often accomplished by capturing, invasive marking and re- 
capturing (DONNELLY et al. 1994, Davis & OVASKA 2001, 
WILLIAMS et al. 2002). However, those invasive markings, 
e.g., toe-clipping in amphibians (PHILLOT et al. 2007) or 
passive integrated transponder tags in fishes (reviewed in 
MussELMAN et al. 2017), may influence the behaviour or 
survival rate of individuals (reviewed in WILSON & MCMa- 
HON 2006). For example, visible implant tags negatively af- 
fected European fire salamander larvae, leading to lower 
body condition indices of marked larvae than newly cap- 
tured larvae (WAGNER et al. 2020a). In recent years, non- 
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invasive methods have gained importance to minimise the 
potential effects on the captured animals. Individual ani- 
mals can be identified by genotyping DNA from faeces or 
hair samples (e.g., BERRY et al. 2014, KUBASIEWICZz et al 
2017), by acoustic surveys such as analysing the distinctive 
howls of Indian wolves (SADHUKHAN et al. 2021) or by pho- 
tographs of individual colouration or patterns, e.g., with 
camera traps (reviewed in CHOO et al. 2020). 

Past photographic mark-recapture studies often includ- 
ed “by-eye” matching of photographs and were performed 
by more than one observer, probably resulting in high er- 
ror rates and variability among observers (e.g., TREILIBS et 
al. 2016). The development of computer-assisted systems 
for photographic mark-recapture studies has led to sub- 
stantial decreases in error rates and between-observer var- 
iability (CRUICKSHANK & SCHMIDT 2017), allows the analy- 
sis of large datasets (BOLGER et al. 2012) and is applicable 
for different species (e.g., JACKSON et al. 2006, VERNES et al. 
2009, DE GASPERIS et al. 2017, RENET et al. 2019). A recent 
study has shown that photographic recognition softwares 
can also be applied to the less suspicious European fire sal- 
amander larvae (FAUL et al. 2022). The non-commercial, 
web-based Wildbook software even combines scientific 
research, crowdsourcing and citizen science and includes 
specialised platforms for a wide range of animals such as 
zebras, sea turtles, whales and giraffes (BERGER-WOLF et al. 
2018). Another recently launched platform, the Amphibi- 
an and Reptile Wildbook (available from www.amphibian- 
reptile.wildbook.org), can currently be used for the pho- 
tographic identification of amphibians such as adult yel- 
low-bellied toads (Bombina variegata) and adult and larval 
European and Near eastern fire salamanders (Salamandra 
salamandra and S. infraimmaculata). 

The European fire salamander occurs across central, 
west and south Europe and is categorised as least con- 
cerned along this range, although the population trend is 
decreasing (KUZMIN et al. 2009). In Germany, which holds 
about 10% of the species range, the status of the Europe- 
an fire salamander has recently been updated to be on the 
pre-warning list (SCHLUPMANN & VEITH 2020). The popu- 
lation decline is probably related to increased human im- 
pacts in forest areas, spring droughts (SCHLUPMANN & 
VEITH 2020) and the introduction of Batrachochytrium 
salamandrivorans (Bsal, MARTEL et al. 2013). Bsal is a path- 
ogenic chytrid fungus that causes the erosive skin disease 
chytridiomycosis in European fire salamanders (MARTEL 
et al. 2013) and has led to population collapses in the Neth- 
erlands, Belgium and parts of Germany (SPITZEN-VAN DER 
SLuys et al. 2016, STEGEN et al. 2017, WAGNER et al. 2019, 
LOTTERs et al. 2020, SCHULZ et al. 2020). Thus, population 
monitoring is crucial to detect potential disease outbreaks, 
prevent the spread of Bsal and maintain stable populations. 

In this study we used the Amphibian and Reptile Wild- 
book to estimate the population size and investigate the in- 
dividual growth of larval European fire salamanders based 
on count data from regular monitorings including photo- 
graphic marking in the Kottenforst, a broadleaf forest near 
Bonn, Germany. Furthermore, we tested for the presence 


180 


of Bsal to confirm its absence since the last Bsal-negative 
reports in the study area (LOTTERs et al. 2018). Using linear 
mixed effect models, we investigated potential effects of the 
larval habitat on population size (counts and estimations), 
mean larval size, individual growth, percentage of injured 
larvae, apparent survival and recapture rates. Although 
the Kottenforst salamander population has been subject to 
many studies investigating differences between the pond 
and stream habitat and potential speciation processes (e.g., 
STEINFARTZ et al. 2007, HENDRIX et al. 2017, SANCHEZ et al. 
2017, SABINO-PINTO et al. 2019, OSWALD et al. 2020), their 
population size is unknown. 


Methods 
Study area and species 


The European fire salamander is a biphasic amphibian 
with aquatic larvae and terrestrial adults, but during the 
reproductive phase, adult females return to water bod- 
ies to deposit their larvae. Usually, females prefer first-or- 
der streams for larval deposition (THIESMEIER 2004), but 
they have also been observed to breed in ephemeral ponds 
and even in caves (e.g., STEINFARTZ et al. 2007, MANEN- 
TI & FICETOLA 2013, CASPERS et al. 2015, MANENTI et al. 
2017). In our study area, the Kottenforst (50°39'38.9” N, 
7°0416.7” E), a forest near Bonn, Germany, females use 
both ponds and streams for reproduction and this behav- 
ioural differentiation corresponds to the differentiation 
into a pond and a stream genetic cluster (HENDRIX et al. 
2017, STEINFARTZ et al. 2007). 

From 2019 to 2022 we monitored European fire sala- 
mander larvae (Salamandra salamandra) and took Bsal 
samples via skin swabs of adult and larval salamanders in 
the Kottenforst. Initially, we chose 5 ponds and a section of 
20 m in each of 5 streams. The exact location of the ponds 
and streams were marked via GPS-coordinates (Fig. 1) and 
the 20m section of each stream was marked at the begin- 
ning, the middle and the end. This enabled us to always 
monitor the same section of the stream. 


The Amphibian and Reptile Wildbook 


The Amphibian and Reptile Wildbook (available for free 
from https://amphibian-reptile.wildbook.org/) is a new- 
ly launched web-based software developed by Wild Me 
(https://www.wildme.org) and is based on computer vision 
algorithm and deep learning to reliably detect and identi- 
fy individual animals from photos. So far, the software is 
trained for adult yellow-bellied toads and adult and larval 
European fire salamanders (SCHULTE et al. 2022). The soft- 
ware allows the upload of large amounts of photographs 
without prior processing. It is also possible to add metadata 
such as GPS coordinates, time and date stamps or the size 
of the animal. Furthermore, it combines citizen science and 
scientific research and can thus provide a huge database and 
valid basis for conservation projects (SCHULTE et al. 2022). 


Fire salamander population analyses via photo-recognition software 


Larval monitoring 


We performed the larval monitoring following a set proto- 
col (see supplementary material) with the time period for a 
certain area of the sample site as the limiting factor for the 
monitoring duration. Each pond was monitored for 30 min- 
utes by one person. Each additional researcher reduced the 
total monitoring time, since the monitoring area per person 
was reduced. Larval density in streams is generally lower 
(WEITERE et al. 2004) and the terrain around streams was 
more difficult to traverse. Therefore, the total monitoring 
time for a 20 m section for one person was set to one hour. 
Except for some occasions in 2019, most of the monitoring 
events were performed by two researchers per sample site. 
In the first year, we performed 10 monitoring events 
of 5 ponds and 5 streams. We visited each pond or stream 
weekly during the first 8 monitoring events from March to 
May, one month later followed by two more monitoring 
events in June and in July. Due to the Corona pandemic, in 
2020 we performed 3 monitoring events (once per month) 
of the same 5 streams and 4 ponds (one pond desiccated in 
2019) in April, May and June. The monitoring period was 
elongated in 2021 to a total number of 15 occasions, aiming 


e pond 
è stream 


to get year-round population estimates of the larval Euro- 
pean fire salamander population. That year, the monitoring 
was performed weekly from March to May, then monthly 
from the end of May to beginning of December. In 2022, 
we started with a monthly monitoring in January and Feb- 
ruary. We then performed 9 weekly monitoring events 
during the field season from March to May and continued 
with monthly monitoring events from end of May until 
September. Due to the low amount of water, we excluded 
one stream (AB) from the monitoring. One pond location 
(TG) had to be replaced by another pond (WT), as it was 
already desiccated in the beginning of the monitoring pe- 
riod. The monitoring events during the four years are visu- 
alised in Fig. 2 (see also Supplementary Table S1). 

During the monitoring of the ponds, we first checked 
visually whether there were larvae visible. Visible larvae 
were caught by gently moving a dip net underneath the 
larvae and then slowly moving it upwards to the water sur- 
face. If no larvae were seen, we systematically went through 
the water with dip nets to capture the European fire sala- 
mander larvae. In ponds, we moved the dip net from the 
centre of the pond towards the littoral zone, because lar- 
vae tended to hide under the leaves and branches or float 


Created with QGIS Version 3.22.5. 
Map base: Humanitarian under CC BY 3.0. 


Figure 1. Map of the study area, the Kottenforst, near the city Bonn, Germany. The inset map of the Kottenforst on the upper right 
shows the ponds (blue circles) and streams (yellow circles) that were monitored regularly for four years. 


181 


Pra Oswa_p et al. 


near the surface in the littoral zone. Once, we could not 
capture any more larvae at one place (i.e., if several strokes 
with the dip net did not return any larvae), we moved on to 
the next part of the pond until the sampling time was over. 
In streams, we again first checked for larvae visually. Once 
detected, we positioned the dip net in an upright position 
downstream of the larvae and moved the dip net towards 
the larvae while simultaneously blocking the way upstream 
with one hand wearing nitrile gloves to prevent touching 
the sensitive skin of larva with our hands. If there was no 
visual detection of larvae, we moved the dip net upstream 
with small up and down movements near the ground, as 
most larvae were sitting on the ground instead of floating 
on the surface. We monitored a section of 1m for approxi- 
mately 3 minutes and moved the dip net very carefully to 
not disturb the environment and organisms living in the 
stream. After about 3 minutes, we moved on to the next 
meter of the stream. We used an interval training App to 
stick to the time periods. For both ponds and streams, we 
put the captured larvae into an opaque 10 | bucket filled 
with water from the respective water body for a maximum 
time of one and a half hours. Afterwards, we transferred 
each larva into a Petri dish with millimetre paper as a scale 
underneath to measure the total body length (from snout 
to tail, in cm, according to Krause et al 2011). Then, we 
transferred it into a small transparent tank with tap water 
(to reduce turbidity of water) to take a photograph from 
the right side of each larva with a digital camera (Sony Cy- 
bershot DSC WX220) or a mobile phone camera (Xiaomi 
Redmi Note 8 Pro, Fairphone 3). We noted specific fea- 
tures such as missing legs or tails. The handling time (size 
measurement and photograph) took maximal three min- 
utes per individual. Afterwards, we collected the larvae in 
an enclosure separate from the bucket to prevent repeated 
sampling. We released all larvae to their capture site im- 
mediately, when the sampling was finished. In addition, we 
randomly sampled larvae from each location each year for 
Bsal (Supplementary Table S2). Bsal was not present at our 
study site (OswaLp et al. unpublished). 


Population analyses 


We uploaded all photographs of the larvae to the Amphib- 
ian and Reptile Wildbook via bulk import clustered for 
each sample site and monitoring event. We added infor- 
mation about the sampling site, date, id, observer and lar- 
val size (as snout-to-tail length). The website uses a specific 
algorithm to detect and match European fire salamander 
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larvae a, if they have been captured in former monitoring 
events. Wildbook does not require any processing of the 
images prior to the detection and matching process. We 
ran the matching process to analyse the number of recap- 
tures and used the Wildbook output (capture history in bi- 
nary format, e.g., 01101, where o = no capture and 1 = cap- 
ture) for population analyses using the Jolly-Seber model 
(POPAN formulation) for open populations as described 
in WAGNER et al. (2020a). This model calculates survival 
probabilities, capture probabilities, the probability of new 
entries (e.g., via birth) and populations sizes. The model 
also allows for emigration by drift and metamorphosis and 
accounts for the prolonged breeding season of amphib- 
ians (WAGNER et al. 2020a). We used the R (R Core Team 
2020) packages dplyr (WICKHAM et al. 2022) and magrittr 
(BACHE & WICKHAM 2022) for data processing and R2u- 
care (GIMENEZ et al. 2018) to test for goodness of fit. Since 
some models based on monthly monitoring data lacked 
goodness of fit for some sample sites, we calculated the ap- 
parent survival rate, the superpopulation size and the esti- 
mated population size including standard errors based on 
weekly monitoring data (March to May) with the package 
RMark (LAAKE 2013). Monthly monitoring data was partly 
used for the boxplots, but not for further analyses or mod- 
els, since we would lose datapoints from multiple weeks 
that were averaged for monthly estimates. 

We calculated the individual daily growth rates (g) 
based on the data of recaptured larvae in 2021 (March to 
December) and 2022 (January to September) by dividing 
the change in size between two successive captures by the 
number of days between the two captures: 

S-S 
p 
where s, = snout-to-tail length at the recapture r, s, = 

snout-to-tail length at previous capture c and t = number 
of days between the two captures. This formula was modi- 
fied from a study on European fire salamanders by Limon- 
GI et al. (2015), who used the same formula but with weight 
data instead of the snout-to-tail length. 

The recapture rates at a given monitoring (r,) were cal- 
culated by dividing the number of recaptured individuals 
by the total number of captured individuals (e.g., JAMES- 
PIRRI & COBB 2000, DEROBA et al. 2005) with the follow- 
ing formula: 


m; 


> 


k on 
k? 
where m, = number or recaptured individuals at a sam- 
pling event k and n, = total number of captured individuals 
at a sampling event k. 


MULLI A LIWALI 
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2021 2022 


Figure 2. Overview of the monitoring events in the years 2019 to 2022. The abbreviations within the horizontal arrow indicate the 
months January to December. Small black arrows indicate a monitoring event. 
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Fire salamander population analyses via photo-recognition software 


Table 1. Linear (mixed) models used for the analysis with the response variable, random effects, fixed effects and their corresponding 
p-values. Significant values are in bold. LME = Linear mixed effects model, LM = Linear model. Due to non-normal distribution some 
variables were transformed. The respective transformation is stated underneath each response variable in parentheses. 


Model Response variable Random effect Fixed effect p-value 

LME1 Number of observed larvae Sample site/year Year p = 0.007 
(ordernorm transformation) 

LME2 Water temperature Sample site/year Habitat type p = 0. 603 
(ordernorm transformation) 

LME3 Mean larval size Sample site Habitat type p = 0.172 

(ordernorm transformation) Session/year Water temperature p = 0.928 

Year p = 90.150 

LME4 Individual daily growth rate ID Habitat type p = 0.907 

(ordernorm transformation) Timespan Year p = 0.382 

Water temperature p = 0.506 

LM5 _ Recapture rate = Habitat p = 0.653 

LME6 Number of estimated larvae Sample site Year p = 0.019 

Session/year Habitat type p = 0.474 

LM7 Estimated superpopulation size - Year p = 0.570 


The percentage of injured larvae (i,) was calculated by 
dividing the number of injured larvae (e.g., missing limbs 
or tailfin, bite marks) by the total number of captured in- 
dividuals: 


b 
i= 
k Nn, ; 
with b, = number of injured larvae at sampling event k 


and n, = total number of captured larvae at sampling 
event k. 

Plots and figures were created with the packages ggplot2 
(WICKHAM 2016), ggpubr (KASSAMBARA 2020), grid 
(MuRRELL 2002) and gridExtra (AUGUIE 2017) in R. 


Statistical models 


We ran linear models with the packages Ime4 (BATEs et 
al. 2015) and ImerTest (KUZNETSOVA et al. 2017) in R (see 
Table 1) to analyse which factors (e.g., habitat type, year) 
affect the mean larval size, the individual growth rate, the 
percentage of injured larvae, the recapture rate, the ap- 
parent survival rate and the estimated number of larvae. 
We also ran a linear model to check whether both habi- 
tats, pond and stream, differ in their temperature regime, 
as amphibians are ectotherm and thus highly susceptible to 
external temperatures. The exact models will be described 
in more detail below. Variables based on population analy- 
ses with the recapture data from Wildbook were calculated 
for the year 2021 and 2022, while the number of captured 
larvae, mean larval size and the percentage of injured lar- 
vae was based on raw data from 2019-2022. For compa- 
rability between years, we only used data from March to 
mid-May. The stream KoB was excluded from the models, 
since it revealed characteristics of both, ponds and streams, 
and could thus not be clearly assigned to any of the two 
habitat types. 


Prior to model selection, we inspected the histogram 
of the variable of interest. We further performed a Sha- 
piro Wilk test and variance test (F test for normal distrib- 
uted data or Fligner test in case of non-normal distribu- 
tion) to check for normality and homogeneity of variances 
(Supplementary Table S3). In case of non-normal data 
distribution, we used the bestNormalize package (PETER- 
SON & CAVANAUGH 2020, PETERSON 2021) to choose the 
most suitable transformation for the data which in most 
cases was done with ordernorm transformation. After- 
wards, we again performed a Shapiro Wilk and variance 
test. For model selection, we set up models with the vari- 
ables of interest including different fixed and random ef- 
fects and chose the model based on the lowest AIC or best 
model performance score as implemented in the perfor- 
mance package (LUDECKE et al. 2021) in R. To validate the 
suitability of the chosen model, we inspected the diagnos- 
tic plots with a focus on the normality of residuals of the 
model with the performance package. If the transformed 
data did not normalise the data, we set up models with the 
non-transformed data using normal and other distribu- 
tions (e.g., Poisson). We investigated the model residuals 
with diagnostic plots and checked for dispersion with the 
DHARMa package (HARTIG 2020) in R. If the diagnostic 
plots and model checks confirmed the underlying assump- 
tions (e.g., normality of residuals), we used the most sup- 
ported model. In case the assumptions were not confirmed 
we performed a simple non-parametric Wilcoxon test (also 
Wilcoxon rank sum test or Mann-Whitney U test) for a 
comparison of means to investigate potential differences 
between larvae from the two habitats (pond and stream). 

After data inspection and model checks, we ran sepa- 
rate linear mixed (effect) models (LM and LME) for the 
dependent variables water temperature, mean larval size, 
individual daily growth rates and recapture rates. For de- 
tailed information on the models see Table 1. The percent- 
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age of injured larvae and apparent survival rates were com- 
pared between habitats with a Wilcoxon test, because the 
model assumptions were not met. 


Results 
Number of captured larvae 


The total number of captured larvae per year ranged from 
406 to 2822 larvae and differed between single monitor- 
ing events within the years, between the years and between 
sample sites (Fig. 3, Table 2, Supplementary Table S1). For 
all sample sites, the number of captured larvae reached 
its maximum between mid-April to mid-May and then 
decreased until August. Due to desiccation, not all water 
bodies were monitored the same number of times. We ex- 
cluded the stream (AB) from the monitoring due to the 
absence of water. One pond (TG) was desiccated after only 
two monitorings and we thus replaced it by another pond 
(WT). 
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Table 2. Number of monitoring events and number of captured 
larvae from 2019 to 2022. The mean number of captured larvae 
(and its standard deviation = SD) is calculated across all sample 
sites per monitoring event and year. 


Year Number of Number of Mean number of 
monitoring events captured larvae captured larvae 
(+SD) 
2019 10 2822 32.67 (+21,64) 
2020 3 406 27,77 (+18.13) 
2021 15 1439 20.36 (+13.72) 
2022 15 2017 26.30 (+32.00) 


Using the dataset including all years and sample sites 
(except for KoB), there was a significant influence of the 
year on the observed number of larvae (LME1, p = 0.007). 
The overall number of larvae decreased from 2019 to 2021, 
but there was a slight increase in 2022 (Supplementary 
Fig. S1). 
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Figure 3. Comparison of the number of larvae per sample site that have been captured during the monitoring events in 2019, 2020, 
2021 and 2022. Shown here are only monitoring events that were conducted during similar time periods in the three years from end 
of March until end of July. Due to the low number of monitoring events in 2020, there are maximum three data points in April, May 
and June per sample site. Interrupted lines or lines that do not reach until August have been desiccated. The pond TG was replaced 


by the pond WT in 2022, since it was completely desiccated. 
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Fire salamander population analyses via photo-recognition software 


Table 3. Range of mean larval sizes (in cm) of larvae captured from ponds and streams for the years 2019 to 2022. The range includes 
the minimum and maximum larval size per year and habitat type. The mean larval size (and its standard deviation = SD) is averaged 
across all sample sites per habitat type and year. 


Year Larval size range Mean larval size (tSD) Larval size range Mean larval size (tSD) 
in ponds (cm) in ponds (cm) in streams (cm) in streams 

2019 2.82-4.76 3.61 (+0.44) 2.75-4.78 3.57 (+0.46) 

2020 3.83-4.11 3.96 (+0.14) 3.23-4.08 3.66 (+0.34) 

2021 2.96-4.40 3.63 (+0.39) 2.93-4.61 3.51 (+0.45) 

2022 3.13-5.01 3.69 (+0.38) 3.00-5.13 3.62 (+0.54) 


Mean larval sizes 


(Table 1). The water temperature did not differ between 


habitats (LME3, p = 0.603, Supplementary Fig. S3). 
Larval sizes varied among monitoring events, sample sites 
and years and ranged from 2.75 to 5.01 cm (Fig. 4, Table 3). 
The mean larval size increased over the course of the moni- 
toring period with larvae captured during later monitoring 
events being larger (Fig. 4, Fig. 5). The habitat type (LME3, 
p = 0.172), water temperature (LME3, p = 0.928) and year 
(LME3, p = 0.150) did not influence the mean larval size 


Individual daily growth rate 


The individual daily growth rate (from first to last cap- 
ture) was calculated for the larvae captured in 2021 and 
2022 and ranged from -12.86 to 15.83% (Table 4). Using 
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Figure 4. The mean larval sizes per sample site during the monitoring events in 2019, 2020, 2021 and 2022. Shown here are only 
monitoring events that were conducted during similar time periods in the three years from end of March until end of July. Due to the 
low number of monitoring events in 2020, there are maximum three data points in April, May and June per sample site. Interrupted 
lines or lines that do not reach until August have been desiccated. The pond TG was replaced by the pond WT in 2022, since it was 
completely desiccated. 
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Table 4. Range of individual daily growth rate (in %) of larvae captured from ponds and streams for the years 2021 and 2022. The 
range describes the minimum and maximum individual daily growth rate. The mean (and its standard deviation = SD) is calculated 


across all sample sites per habitat type and year. 


Year Range of growth rates (%) Mean growth rate (%) Range of growth rates (%) Mean growth rate (%) 

in pond larvae in pond larvae (+SD) in stream larvae in stream larvae (+SD) 
2021 -0.18-7.14 1.52 (£1.21) - 2.14-4.29 1.01 (+1.23) 
2022 -4,00-11.67 1.28 (+1.80) -12.86-15.83 1.19 (+3.41) 


the full dataset, there was no effect of the larval habitat 
type (LME4, p = 0.907), the year (LME4, p = 0.382) or 
the water temperature (LME4, p = 0.506) on daily growth 
(Fig. 5). 


Percentage of injured larvae 


The percentage of injured larvae (limbs or tailfin miss- 
ing, bite marks) across all years ranged from o to 100% in 
pond larvae (mean 11.84%) and o to 37.5% in stream larvae 
(mean 9.80%). There was no significant difference between 
the two habitat types (Wilcoxon, p = 0.250, Fig. 6), but a 
significant difference between the years (Kruskal-Wallis, 
p < 0.001). When we split the dataset into the single years 


Daily growth rate (cm) 


Pond 


(Fig. 7, Table 5), there was still no significant difference in 
the percentage of injured larvae from the two habitat types 
in 2019 (Wilcoxon, p = 0.069) and 2020 (LMs, p = 0.653), 
but there was a significance effect of the habitat type in the 
two following years. In 2021, there was a higher percentage 
of injured larvae in streams (Wilcoxon, p = 0.016). In 2022, 
ponds exhibited a higher number of injured larvae (Wil- 
coxon, p = 0.037). 


Recapture rates 
Recapture rates were based on the photographic recap- 


tures during the monitoring events from March to May 
in 2021 and 2022. The recapture rates varied between sam- 


Stream 


Habitat type 


Figure 5. Individual daily growth rate per habitat type across all sample sites (excluding KoB) based on the monitoring events in 2021 
and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the median (black horizontal 
line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 times the interquartile 
range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as outliers. The annotation 
above the box plot represents possible differences of individual daily growth rates between fire salamander larvae from the two habitat 
types based on p-values: ns = non-significant. 
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Figure 6. Percentage of injured larvae per habitat type across all sample sites (excluding KoB) during the monitoring events in 2019, 
2020, 2021 and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the median (black 
horizontal line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 times the 
interquartile range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as outliers. 
The annotation above the box plot represents possible differences of the percentage of injured larvae between the two habitat types 
based on p-values; ns = non-significant. 
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Figure 7. Percentage of injured larvae per habitat type across all sample sites (excluding KoB) per year during the monitoring events 
in 2019, 2020, 2021 and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the median 
(black horizontal line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 times 
the interquartile range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as outli- 
ers. The annotation above the box plot represents possible differences in the recapture rate between the two habitat types based on 
p-values: ns = non-significant, * = p < .05 
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Table 5. Range of injured larvae (in %) from ponds and streams in the years 2019 to 2022. The range describes the minimum and 
maximum percentage of injured larvae in ponds and streams. The mean (and its standard deviation = SD) is calculated across all 
sample sites per habitat type and year. 


Year Range of the number of Mean number of injured larvae Range of the number of | Mean number of injured larvae 
injured larvae in ponds (%) (%) in ponds (+SD) injured larvae in streams (%) (%) in streams (+SD) 

2019 0-20.59 3.52 (+5.43) 0-27.27 5.68 (+6.31) 

2020 0-20 7,43 (+9.55) 0-30 9,11 (410.73) 

2021 0-23.08 7,18 (+7,48) 0-37.50 13.70 (£10.21) 

2022 0-100 25.51 (+24.87) 0-30 10.71 (+8.71) 


Table 6. Range of recapture rates (in %) for larvae from ponds and streams in the years 2021 and 2022. The range describes the mini- 
mum and maximum recapture rates (in %). The mean recapture rate (and its standard deviation = SD) is calculated across all sample 
sites per habitat type and year. 


Year Range of recapture rates Mean recapture rates (%) Range of recapture rates Mean recapture rates (%) 


in ponds (%) in ponds (+SD) in streams (%) in streams (+SD) 
2021 0-50 11.65 (+14.30) 0-100 8.01 (+9.06) 
2022 0-50 9.51 (411.87) 0-58.78 10.47 (412.14) 


ple sites and monitoring events (Supplementary Table S4, in streams. In 2022, the recapture rates for monitoring 
S5). In 2021, the highest recapture rate of 100% (ina pond) events with sample sizes higher than ten were up to 37.70% 
was found during a monitoring event witha very low sam- in ponds and 58.70% in streams (Table 6). The recapture 
ple size of only one individual. During monitoring events rate (including both years) did not differ between habi- 


with sample sizes higher than ten, recapture rates were tats (Wilcoxon, p = 0.780, Fig. 8) or years (Wilcoxon, p = 
still high and reached up to 24.14% in ponds and 34.04% 0.898). 
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Figure 8. Recapture rate per habitat type across all sample sites (excluding KoB) during the weekly monitoring events from March 
to May in 2021 and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the median 
(black horizontal line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 times 
the interquartile range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as outli- 
ers. The annotation above the box plot represents possible differences in the recapture rate between the two habitat types based on 
p-values: ns = non-significant. 
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Table 7. Range of the apparent survival (in %) of larvae from ponds and streams in the years 2021 and 2022. The range describes the 
minimum and maximum survival rates (in %). The mean (and its standard deviation = SD) is calculated across all sample sites per 


habitat type and year. 


Year Range of the apparent Mean apparent survival (%) Range of the apparent Mean apparent survival (%) 
survival in ponds (%) in ponds (+SD) survival in streams (%) in streams (+SD) 

2021 48.81-93.67 76.58 (+17.45) 79.01-100 89.24 (+8.79) 

2022 0-100 57.89 (+44.87) 0-100 64.77 (+43.06 


Survival rates 


The apparent survival rates were based on the photo- 
graphic recapture data during the weekly monitorings 
from March to May in 2021 and 2022 and ranged from 
o to 100% (Table 7). There was no influence of the year 
on the apparent survival (Wilcoxon, p = 0.502), but the 
apparent survival differed significantly between habitat 
types (Wilcoxon, p = 0.034, Fig. 9). 


Population size estimation 


The estimated population size was based on the photo- 
graphic recapture data from 2021 and 2022 and was high- 
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ly variable between sample sites and monitoring events 
(Fig. 10). Weekly and monthly estimates were similar in 
2021, but deviated in 2022 (Supplementary Fig. S3, S4). 
In the following, we will focus on weekly estimates, since 
the time periods between weeks were more regular and 
thus more suitable for population analyses. The estimat- 
ed number of larvae in 2021 ranged from 43.75 to 1330.30 
individuals in ponds and from 15.34 to 561.00 individu- 
als in streams (Supplementary Table S6). In 2022, the es- 
timated larval population size ranged from o to 972.92 
individuals in ponds and from o to 647.97 individuals 
in streams (Table 8, Supplementary Table S7). The year 
(LME6, p = 0.019), but not the habitat type (LME6, p 
=0.474) had a significant influence on the number of es- 
timated larvae. 


Stream 


Habitat type 


Figure 9. Apparent survival rate per habitat type across all sample sites (excluding KoB) during the weekly monitoring events from 
April to May in 2021 and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the median 
(black horizontal line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 times 
the interquartile range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as outli- 
ers. The annotation above the box plot represents possible differences of the survival rate in the two habitat types based on p-values: 


ns = non-significant, * = p < .05. 
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Table 8. Range of the estimated larval population size. The range describes the minimum and maximum estimated number of larvae 
per habitat type and year (including the standard deviation = SD). 


Year Range of the estimated number of larvae in ponds (+SD) Range of the estimated number of larvae in streams (+SD) 
2021 43.75 (£15.92) — 1330.30 (+690.00) 15.34 (+1.16) - 561.00 (+174.90) 
2022 0 (+0) - 972.92 (+35.93) 0 (+0) - 647.97 (+60.90) 


The estimated size of the superpopulation in 2021 
ranged from 172.94 (+46.32) to 4495.61 (+2084.50) individ- 
uals in ponds and from 383.32 (+137.36) to 1013.93 (+239.18) 
in streams. In 2022, the estimated superpopulation size 
ranged from 67.81 (411.83) to 1264.90 (+13.05) in ponds 
and from 480.76 (+0.00) to 16236.24 (+2622.44) in streams 
(Supplementary Table S8). There was no significant differ- 
ence in the estimated superpopulation size between the 
years 2021 and 2022 (LMz, p = 0.570, Fig. 11). 


Discussion 


In our study we investigated mean larval size, recapture 
rate, apparent survival, estimated larval population size, 
individual growth and the percentage of injured larvae of 
a European fire salamander population in Germany. We 
found high variation in those traits within and between 
sample sites and monitoring events. There was an effect of 
the year but not of the habitat type on the observed number 
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Figure 10. Larval population size estimates across all sample sites in the Kottenforst based on weekly monitoring events in 2021 and 
2022. The upper row represents the estimated number of larvae for the streams, the bottom row represents estimated population sizes 
for the ponds. Stream KoB is in parentheses, since it has characteristics of both, pond and stream. In many cases, the monthly estima- 
tion of population size after July was not possible, since there were no more larvae found. Sample sites KoB, TG and SG desiccated 


after April and May, respectively. 
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of larvae, the percentage of injured larvae and the estimat- 
ed population size. There was no difference in superpopu- 
lation size. When split into subsets for each year, the per- 
centage of injured larvae differed between habitat types in 
2021 (higher in streams) and 2022 (higher in ponds). There 
was no effect of the year or habitat type on recapture rates. 
Recapture rates, mean larval sizes and water temperature 
were independent of the habitat and year. The mean larval 
size was not influenced by the water temperature. The ap- 
parent survival rate and the daily growth rate differed sig- 
nificantly between habitats, although this effect was prob- 
ably driven by year-specific differences. We found no cases 
of Bsal in our study area. 

So far, only few studies have investigated European fire 
salamander population sizes (but see SEGEV et al. 2010, 
ALVAREZ et al. 2015, SINAI et al. 2020). According to two 
previous studies, population sizes were higher in perma- 
nent water bodies (SEGEV et al. 2010, SINAI et al. 2020) 
indicating permanent water bodies provide a more suit- 
able habitat than temporary water bodies. However, in our 
study there was no significant difference between tem- 
porary ponds and permanent streams regarding the esti- 
mated population size. This discrepancy might result from 
different monitoring approaches. SEGEV et al. (2010) and 
SINAI et al. (2020) monitored adults in proximity to the 
water bodies, while we monitored larvae in the respective 
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habitats. Nevertheless, we found a higher apparent surviv- 
al in streams than in ponds, indicating that streams, i.e., 
here the permanent water body, might indeed be advan- 
tageous. 

We found no significant difference in the individual 
daily growth rate across the overall dataset with regard to 
the habitat type. These findings are in line with a previ- 
ous study that did not find growth difference between lar- 
vae from the two habitat types (SABINO-PINTO et al. 2019). 
However, during a reciprocal transfer experiment that we 
conducted in 2019, we found higher growth rates before 
metamorphosis in larvae that were transferred to ponds 
(OswALD & CASPERS, unpublished). Other studies found 
that larval growth was affected by temperature, desicca- 
tion events and food regime (e.g., ALCOBENDAS et al. 2004, 
Krauss et al. 2011). While low food abundance led to de- 
creased growth rates in fire salamander larvae (e.g., ZAHN 
2007, LIMONGI et al. 2015), desiccation risk accelerated lar- 
val development at the expense of reduced growth (RicH- 
TER-Borx et al. 2011, SZEKELY et al. 2017). Although we 
did not find temperature differences between ponds and 
streams in our study, both habitats were different with re- 
gard to food abundance and desiccation risk (P. OswaLp, 
unpubl. obs.). As we did not find overall differences in 
growth, larvae seem to be able to compensate for disad- 
vantageous conditions in the larval habitat. 


2022 
Year 


Figure 11. Larval superpopulation size per year across all sample sites (excluding KoB) based on the weekly monitoring data from 
March to May in 2021 and 2022. The box plots are representing the upper and lower quartile (upper/lower edge of the box), the me- 
dian (black horizontal line within the box), and the maximum and minimum values (end of upper/lower whisker) lying within 1.5 
times the interquartile range. Data points that are > 1.5 times the interquartile range (upper quartile-lower quartile) are defined as 
outliers. The annotation above the box plot represents possible differences in the estimated number of larvae between the two years 


based on p-values: ns = non-significant. 
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Previous studies often considered streams as the more 
suitable habitat for European fire salamander larvae (e.g., 
THIESMEIER 2004, REINHARDT 2014). In concordance with 
this, we found higher survival rates in streams. We have 
to consider that our approach calculates apparent survival, 
i.e., we cannot distinguish between emigration and death 
(LEBRETON et al. 1992). Apparent survival might thus be 
misleading and does not provide true survival estimates 
(often underestimated), since emigrated individuals might 
still survive outside of the investigated population (Gr1- 
ROY et al. 2012). As larval drift has a large influence on sur- 
vival and catastrophic events such as strong precipitation 
decreases redetection probability (REINHARDT et al. 2018), 
downstream drift is likely to result in underestimated sur- 
vival estimates in streams. In our study, the apparent sur- 
vival for stream larvae was quite high indicating that either 
survival estimates were less underestimated than expect- 
ed or survival was even higher. Since population estimates 
did not differ significantly between both habitats (but see 
comment above) and the apparent survival was quite high 
in both habitats, ponds and streams might provide equally 
suitable habitats. 

Population estimates such as survival or population siz- 
es are highly dependent on the used methods. For example, 
encounter surveys might best reflect relative abundances, 
while mark-recapture studies facilitate the estimation of 
population parameters such as survival and recapture rates 
(Bower et al. 2014). According to SCHMIDT (2004), sim- 
ple count data neglects potential differences in detection 
probabilities and does not reflect true population sizes. 
This was confirmed by a recent study comparing four dif- 
ferent monitoring methods which found that count meth- 
ods and removal sampling underestimated population 
sizes (WAGNER et al. 2020a). Furthermore, the estimated 
population size is dependent on the breeding system. In 
case of amphibians, which include “explosive” breeders and 
prolonged breeders (WELLS 1977), the superpopulation ap- 
proach might be a useful tool to estimate population sizes 
in open populations where individuals occur asynchro- 
nously during the breeding season (PELLET et al. 2011). 

In accordance with WAGNER et al. (2020b), our study 
emphasizes the need for regular population monitoring, as 
there was high variation between and within years and be- 
tween different locations. Thus, sporadic monitoring does 
not always represent the actual population size and sta- 
tus. Especially, in the context of emerging infectious dis- 
eases such as Bsal, irregular monitoring events might fail 
the early detection of population declines. The choice of 
the appropriate monitoring method, depending on the ob- 
jectives of the study and available resources, is crucial to 
get appropriate population estimates and detect potential 
declines (STORFER 2003). Alternatively, to our mark-cap- 
ture-recapture approach, larval removal sampling, which 
does not require individual identification and can easily be 
conducted in the field, also reflects a less time consuming 
and inexpensive approach to detect population changes 
(SCHMIDT et al. 2015, WAGNER et al. 2020b). Although we 
did not find Bsal in our population, yet, other populations 
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have suffered from harsh declines due to this pathogenic 
fungus (e.g., SPITZEN-VAN DER SLUYs et al. 2013, STEGEN et 
al. 2017, DALBECK et al. 2018). A recent study even found a 
new disease that threatens European fire salamander pop- 
ulations (MANENTI et al. 2022). 

Thus, due to the multiple drivers of amphibian pop- 
ulation declines that are often spatially different, local 
monitoring programs are key to amphibian conservation 
(GRANT et al. 2016). These programs often include capture- 
mark-recapture studies that come with a trade-off between 
the need for demographic estimates, logistic issues and po- 
tential harm to individuals during marking (HADDAD et 
al. 2008). For amphibians with distinctive colour patterns, 
there are several photographic identification softwares that 
vary in their efficacy to recognise and identify individuals 
(MaTTHE et al. 2017). Monitoring efforts are often unco- 
ordinated, unequally distributed among taxa and there is 
a general need for global open databases (Moussy et al. 
2022). The Amphibian and Reptile Wildbook used in this 
study has the potential to become such a database for am- 
phibians and reptiles, as it is easy to use, combines citizen 
science and scientific projects and even works for amphib- 
ian larvae. Larval monitoring has many advantages (easier 
to find, defined area), especially in species where terrestrial 
adults are difficult to survey, e.g., because they are cryptic 
or found in low densities (SKELLY & RICHARDSON 2009). 
The Amphibian and Reptile Wildbook is able to detect and 
identify both, larval and adult individuals of the European 
fire salamander (plus yellow-bellied toads and Near East- 
ern fire salamander). It is an open-source software and free 
to use, the software performs similarly or better than com- 
parable photo-identification programs (L. SCHULTE et al. 
under study) and is robust against varying picture quali- 
ties (SCHULTE et al. 2022). It identifies 99.6% of larval re- 
captures from photos and proved to be a reliable tool in 
European fire salamander larvae research (L. SCHULTE et 
al., under study). Thus, the Amphibian and Reptile Wild- 
book might pave the way for future amphibian conserva- 
tion projects based on both adult and larval monitoring in- 
cluding the opportunity of an open database and interna- 
tional collaborations. 


Conclusion 


We found quite high estimates that varied within and be- 
tween sample sites and fluctuated over time. The popula- 
tion estimates decreased significantly from 2021 to 2022, 
but did not seem to differ between pond and stream hab- 
itats. However, we found a higher apparent survival rate 
in streams, compared to ponds. ‘These results indicate that 
ponds and streams might be equally used as larval habitats, 
but streams might still be beneficial. It provides useful in- 
formation for conservation projects such as the establish- 
ment of new breeding ponds or connectivity measures. Al- 
though the chytrid fungus Bsal is spreading, it has not yet 
reached the Kottenforst. However, the population should 
be monitored regularly to detect possible outbreaks. Fi- 
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nally, the larval monitoring and the new Amphibian and 
Reptile Wildbook software have proven to be efficient tools 
for population analyses of European fire salamanders. This 
method can probably be applied to other monitoring pro- 
grams of amphibian adults and larvae, if the animals exhib- 
it distinctive individual patterns and might be a first step 
towards a global database for amphibian and reptile moni- 
toring data that supports (global) conservation efforts. 


Data availability 


Raw data and code are available at the online repositories GitHub 
and Open Science Framework (https://doi.org/10.17605/OSEIO/ 
7BKQD). The monitoring protocol for the European fire salaman- 
der larvae and additional tables can be found in the Supplemen- 
tary material. 
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Supplementary data 


The following data are available online: 


Supplementary Table S1. Overview of the number of larvae cap- 
tured during the monitoring events from 2019 to 2021. 


Supplementary Table S2. Overview of the Bsal samples from 2019 
to 2022. 


Supplementary Table S3. Results of the tests for normal distribu- 
tion and variance of homogeneity for all investigated variables. 


Supplementary Table S4. Number of recaptures among captured 
animals and recapture rates per sample site based on the moni- 
toring events in 2021. 


Supplementary Table S5. Number of recaptures among captured 
animals and recapture rates per sample site based on the moni- 
toring events in 2022. 


Supplementary Table S6. Estimated larval population size (and 
standard error) per sample site based the weekly monitoring 
events in 2021. 


Supplementary Table S7. Estimated larval population size per 
sample site based the weekly monitoring events in 2022. 


Supplementary Table S8. Estimated larval superpopulation size 
per sample site based the weekly monitoring events in 2021 and 
2022. 


Supplementary Figure S1. Observed number of larvae across all 
sample sites during the monitoring events in 2019, 2020 and 2021 
and 2022. 


Supplementary Figure S2. Water temperature per habitat type 
across all sample sites during the monitoring events in 2019, 2020 
and 2021 and 2022. 


Supplementary Figure $3. Larval population size estimates across 
all sample sites in the Kottenforst based on weekly and monthly 
monitoring events in 2021. 


Supplementary Figure $4. Larval population size estimates across 
all sample sites in the Kottenforst based on weekly and monthly 
monitoring events in 2022. 
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